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^— ■») , We consider the transport spectroscopy of a quantum dot with an even number of electrons at 

' finite bias voltage within the Coulomb blockade diamond. We calculate the tunneling current due 

\ to the elastic and inelastic co-tunneling processes associated with the spin-singlet ground state and 

the spin-triplet first excited state. We find a step in the differential conductance at a bias voltage 
' equal to the excitation energy with a peak at the step edge. This may explain the recently observed 

sharp features in finite bias spectroscopy of semi-conductor quantum dots and carbon nanotubes. 
Two limiting cases are considered: (i) for a low excitation energy, the excited state can decay only 
' by inelastic co-tunneling due to Coulomb blockade (ii) for a higher excitation energy the excited 

CO . state can decay by sequential tunneling. We consider two spin-degenerate orbitals that are active 

in the transport. The nonequilibrium state of the dot is described using master equations taking 
spin degrees of freedom into account. The transition rates are calculated up to second order in 
^ ■ perturbation theory. To calculate the current we derive a closed set of master equations for the 

' spin-averaged occupations of the transport states. 
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I. INTRODUCTION 



\ Transport spectroscopy provides information on the excited states of quantum dots.0 Outside the Coulomb block- 
aded regionS|-.( "diamonds"), peaks in the differential conductance can be observed, which linearly shift with the 
gate voltage.0 Such a peak indicates that an excited state has entered the transport voltage window. Recently, De 
' O ■ Franceschi et al. observed new features in the spectroscopy data of a vertical semiconductor quantum dot coupled 
^ ' relatively strong to the leads B Inside several Coulomb diamonds with an even number of electrons they observed a 
step in the differential conductance with a peak at the step edge for a quantum dot with six electrons. In contrast 
to the peaks outside the diamond, this feature depends very weakly on the gate-voltage. This was interpreted as 
being due to inelastic co-tunneling through an excited state of the dot. It was pointed out that the sharpness of the 
step is determined either the temperature or by the broadening of the excited level, both of which are small energy 
scales compared to the widths of the first-order tunneling peaks. This provides an new spectroscopic tool for studying 
excited levels in quantum dots. Similar features were observed in the transport spectroscopy of carbon nanotubes 
obtained by Nygard et. al..u Here also, the gate voltage independent features at finite bias voltage were clearly seen 
in consecutive Coulomb diamonds with an even number of electrons. The observed peaks at the step edges were 
attributed to Kondo-physics. In both experiments a Kondo effect was observed at zero bias voltage in a magnetic 
field. The occurence of Kondo geaks in the differential conductance, although of different origin, indicates a relatively 
strong coupHng to the leads .013 

Motivated by these experiments, we have calculated the co-tunneling current through an excited level of a quantum 
dot using a master equation approach. We investigate the possibility that the observed features arise from changes 
Ch ■ in the nonequilibrium occupations of the states in the dot. Such changes are caused by voltage-dependence of 
, co-tunneling processes. We find that a peak structure in the differential conductance can develop at the onset of 
^ inelastic co-tunneling. We consider the possible regimes defined by the energy of the excited level relative to the 
Q addition energy. Two possible cases are distinguished^ (i) for low excitation energy, the excited state can decay 
O only by inelastic co-tunneling due to Coulomb blockade; (ii) for higher excitation energy the excited state can decay 
by sequential tunneling. We take the spin degrees of freedom into account and derive master equations for the 
nonequilibrium occupations of the |Celevant states, which are averaged over the spin projections. The case (i) has 
recently been considered by M. EtoS without taking the spin degrees of freedom into account. 

The plan of the paper is as follows. In Sec. || we introduce a model for the transport spectroscopy of a quantum 



dot. The low bias elastic co-tunneling current is discussed in Sec. [II and in Sec. iV we consider the different regions 
of the Coulomb diamond at higher bias where inelastic co-tunneling dominates the transport. Our conclusions are 
presented in Sec. 0. 
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II. TRANSPORT SPECTROSCOPY 



We consider a quantum dot connected to two leads L and R and one gate electrode G. The dot is modeled by 
the Hamiltonian H = Hgp + i?ci + ^^ox- The discrete single particle states are labeled by d — 1,2,... and are spin 
degenerate, s =1,1^ 



ds 

where n^s = O'dgO.ds is the number operator of leV|el ds. We describe the electron-electron interactions in the classical 
electron liquid picture (Coulomb blockade model)EJ. The electrostatic energy depends on the total number of electrons 
N = J2ds ''^ds in the dot and on the voltages on the electrodes, which couple to the island through the capacitances 
Cl, Cij, Cg (total capacitance C — J2k=L rg^^)'- 

h,i = \n{n-i)u+[\u- 

\ k=L,R,G 

Here U ~ 1/C is the charging energy (we use units e = k = % = 1). We also include an exchange term that lowers 
the energy of states with parallel spins: 

^ox — — J ^ ndsnd's 

dd's 

The leads are modeled as noninteracting quasi-particle reservoirs with spin degenerate states. For the left lead with 
states / we write Hl = X^is ^i^is where nis — aj^ais and Hr is defined similarly for the right lead (states labeled by 
r). The leads are assumed to be in equilibrium at chemical potential fJ-L = fJ-R + V and /i/j, respectively. They are 
coupled to the dot by tunnel junctions. The tunneling Hamiltonian of the left junction is 



Htl - Y {^TLd + HtlcL^ , HrLd - Y ^idal 

ds Is 



ads 



and Htr. is defined similarly for the right junction {L,l — > R,r). The temperature T is assumed to be small: 
T ^ r, y, 6, where F is the typical broadening of an energy level in the dot due to the coupling to the leads and 5 is 
the relevant level spacing introduced below. 

We consider here the case where the number of electrons on the dot N is even at zero bias voltage. Let us first 
consider the conditions for sequential tunneling transport through the dot. For transitions between the N and iV ± 1 
electron ground states only two orbitals need to be considered. Let |0) denote the N ~ 2 electron ground state and 
d and d' label the and N/2 + T** orbital with level spacing 5 = 8^/2+1 — £n/2- The singlet ground state then 

reads |-/Vo,o) = oj^icj^jO). Here the subscripts in \Ns,s~) indicate the spin S on the dot and its projection Sz in some 
arbitrary z-direction. The doublet ground states with iV ± 1 electrons are denoted by 

l(^-l)i/2,+i/2> = «dTlO)' l(^ + l)i/2,+i/2> =«d'T4T«dilO>' (la) 

l(^-l)l/2,-l/2>=«di|0>, l(A^+l)l/2,-l/2>=4'i«dT4i|0>, (lb) 

The differences between the ground state energies Ej\[ are 

_ ( l\ ^ Ck 

Enn-1 — Epj — En~i = £n/2 + (^^2/^^ ^ ~C^^' 

^ ^ k=LRG 



En+in = Ej^+i — Epf = 6^/2+1 + (n + 2 ) ^ ^ 



k=LRG 



In order to have sequential tunneling transport it must be possible to inject an electron from the left lead and emit 
one to the right lead in a subsequent event: 

^J■L > Enn-1 > ^J■R (3a) 
fJ-L > En+in > fJ-R (3b) 



2 



These transport conditions translate into a stability region in the plane of bias and gate voltages, which is shown in 
Fig. 1^. Let us express the addition energies as a function of the position {V, Vq) in this diagram. In the limit of zero 
bias voltage (i.e. V = in Fig. Q), sequential tunneling transport is possible only at two isolated points: Vq — Vq^^^ 
[Eq. ( ^ ) is satisfied] and Vg = Vq^^'^ [Eq. ( pb| ) is satisfied]. At these degeneracy points the number of electrons 
on the dot can change by one. At the intermediate gate voltages Vq'^~^ < Vq < Vq'^^'^ the number of electrons 
on the dot is stable and equal to N. Sweeping the gate voltage one thus observes the Coulomb blockade oscillations 
with peak spacing 

yN + lN _ yNN-1 ^ {S+U). 

This involves the level spacing 6 because the number of electrons on the dot is even. For a finite, fixed bias voltage 
each of these degeneracy points develops into a gate-voltage transport window, cf. Fig. |l|. Equivalently, the transport 



conditions (pa|) and (3b) require that the bias voltage exceeds a gate- voltage dependent threshold to allow the number 



of electrons on the dot to fluctuate. The threshold lines are (depicted in Fig. 



V^-' ^ i^o - V^"^-') , (4a) 



V^^-' ^ +g {Va - V^^-') , (4b) 



^ {V^-^^^ - Va) , (4c) 



^fl = [Vg -Vg) ■ (4d) 



The region where the N electron ground state is stable, the A^th Coulomb diamond, is thus delimited by 

yNN-l ^ yN+lN ^y^^ ^y ^ yNN-1 ^ yN+lN g^^^ ^^J^^g^^ ^^A^-l ^ y^ ^ yN+lN^ r^^^ 

mal bias voltage that can be applied without lifting the Coulomb blockade is 

Vax =—[VG V(.)=d + U. (5) 

This is achieved at the intersection of the threshold lines (4b) and (^) (the top of the diamond) at gate voltage 
max = V^^-^ + {Cl/C) V^f - V^^^^^). Similarly, V = -V^^^ at the intersection of hues (|^ and (gdj) at 



Vg = y^"*"^^ — {Cl/C) {Vq^^^ — V^^ ^). The addition energies relative to the electro-chemical potentials can now 
be expressed in terms of distance of the bias voltage to the (gate- voltage dependent) thresholds value 

Enn-i =l^R-^ {y^^^ (^g) -y)< Mfl, (6a) 

En+in = + (l - iyL^'"" i^G) -V)> ML. (6b) 

Let us now consider the transport conditions for the A^ electron excited state. We assume that the lowest lying 
excited state is a spin triplem, 

= 4't4tIo), 

with energy Ejqi = A -I- Eiq. The prime indicates the excited state and the excitation energy A = 5 — J is just the 
spacing between the two active levels d and d' reduced by the exchange energy. We concentrate on the change in the 
current due to this excited state. The effect of excited states with higher energy (such as the spin singlet at En + 5) is 
disregarded. In order to have transport by sequential tunneling, the transitions between the A^ -I- 1 and A^ — 1 electron 
ground states and the A' electron excited state must be possible: 

Ml > En'n-1 > M-R., (7a) 
Ml > En+in' > fJ-R- (7b) 
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where 



En'n-1 — Enn-1 + A, i^AT+iTV' — E'at+itv — A (8) 

The regions in the stabihty diagram where the conditions (|7a| ) and (|3b|) are satisfied, are dehmited by another four 
thresholds fines for the bias voltage: 

Vf^~^ iVa) = Vi^^-i (Va + £a) , (9a) 
Vf''+' (Vg) = C^+i (vg ~ ^a) , (9b) 

respectively, where a — L,R. These are simply the regions in Fig. |l| where ( [7a| ) and (^) are satisfied, shifted in gate 
voltage by (C/Cg)A towards the center of the Coulomb diamond. Outside the Coulomb diamond the sequential 
tunneling transitions to and from the N electron excited state gives rise to an increase in the current. By measuring 
the differential conductance a spectroscopy of excited states can be performed.tl Within the Coulomb diamond the 
N electron ground state is stable with respect to sequential tunneling. If the tunnel coupling to the leads is weak, 
the transport is suppressed within the Coulomb diamond. As the tunnel coupling is enhanced transport becomes 
possible as a result of higher order tunneling processes. At low bias voltage, < A, only elastic co-tunneling is 
allowed by energy conservation. Charge is transported through the dot without exciting it. At higher bias voltage, 
V > A, inelastic co-tunneling processes excite the dot and give rise to an additioaal current. Since these processes 
transfer energy between the dot and the leads, the co-tunneling is called inelastic.cl The dot switches between the N 
electron ground and excited state and in certain cases the TV ± 1 electron ground states can be occupied, too. These 
two regimes are considered separately in the next two sections. We point out that the model introduced here can be 
applied to a limited range of bias and gate voltages in an experimental situation where the capacitances and single 
particle energies provide a sufficiently accurate parametrization of the many-body spectrum of the dot. 



III. ELASTIC CO-TUNNELING 



At low bias voltage < A, a tunneling process of second order in the transition amplitude cannot populate the 
excited state. However, in this order of tunneling, transport is possible via elastic co-tunneling without exciting the 
dot. This is the dominant contribution to the current at sufficiently low temperatures and voltages. We dcunot 
consider here the integer-spin Kondo physics o»,|a-small energy scale Tr-, which has recently been observed: we 
assume T,V '3> Tk- It was pointed out in Refs. tallll that elastic co-tunneling is the only process limiting the operation 
of a quantum dot as a single electron transistor when the level quantization 5 is comparable with the charging energy 
U. 

In a co-tunneling process two tunneling events occur in a very short time interval. For instance, an electron is first 
injected into the dot from the left lead and then emitted to the right lead. On this small time scale the two tunneling 
events are coherent. The energy of the total system (dot and leads) is sufficiently uncertain for the process not to 
be completely blocked. However, the transition amplitude for this proc|3ss is inversely proportional to the energy 
difference between the initial and intermediate state of the total system.ll3 The amplitude for transferring an electron 
in the opposite sequence (first emit an electron, then inject one) via another intermediate state must also be added if 
the final state of the system is the same. From the total second order amplitude (or matrix element) one obtains the 
tunneling rate by statistically averaging over all initial and final states of the leads Thus a current can pass through 
the dot although the iV — 1 and iV -I- 1 electron states are only virtually accessible. 

Let us describe the calculation of the tunneling rate for this simple case in order to skip the details of similar 
calculations below. Consider the initial state of the system where two electrons are in level I of the left leads and the 
dot in the ground state, denoted by \NQ^Q)a\^a\^\) where |) denotes the state of the other electrons in the noninteracting 

leads. Applying the tunneling operators H^j^^ and Hj^^d in the two possible sequences, one obtains the final state 
|A^o,o)-F'|) where the operator F is given in Table |. The amplitude for this process is written as v^Mjv^tv where the 
factor \/2 is due to the spin degeneracy of the levels in the left lead. The explicit expression for M^^n is given in 
Table ||. Introducing the tunneling rate 

Tn^n = 27r ^ iM^^jvl' / (£/ - Mi) [1 - / i^r - Mk)] S (ez - Sr) , (10) 

l,r 

we write for the elastic co-tunneling current 
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2r 



(11) 



The sum over all possible states of the leads with the Fermi distribution / (e) = 1/ (l + e"^/^) can be converted into 
an integration and gives 



^ [^dd^ddT {Enn-i,V) + T^'d'^d'd'T (En+in, V) 
+ ^dd'^d'd'^' {Enn-i,En+in, V)] 



(12) 



d') for sequential tunneling into level d 



The same applies to the quantity F^^, = 



Here we introduced the rates F^^ — 27r^Ji;d| S {e — ei) and F^,^, (d - 
and d' and we assumed that these depend only weakly on the energy e. 

^ttJ^i tid^id'S (e — £i) in the last term on the right-hand side of Eq. (12), which takes into account the interference 
between the two possible sequences for the tunneling events ("paths" in energy space); T^^,, F^^,F^^, are defined 
similarly {L,l — + R,r). All these quantities are treated as parameters here. In the expression for the rate (^2|) the 
explicit (V) and implicit dependence [through the addition energies; cf. Eqs. (^] on the voltages is incorporated in 
the functions (explicitly evaluated in the limit T —f 0) 



r{E,V) 



t'{E, E',V) 



de 



/ (£ - ml) [1 - / (£ - A*fl)] 



V 



{s-EY 



it,L - E) ifiR - E) 

/(£-AiL)[l-/(e-Mfl)] 



de- 



\e~E\ 



E' 



(13a) 



(13b) 



1 



\E' -E\ 



( 


In^^ 


-E 


+ 


In^^ 


-E' 


) 






-E 






- E' 





e(F), 



Here Q is the unit step function and liniE'^E t' {E, E' ,V) — t{E,V). The voltage dependence will be important 
below. The rate (^2|) diverges near the edges of the Coulomb diamond where sequential tunneling becomes dominant 
[cf. condition (0)]. Here the second order perturbation theory for the amplitudes breaks down. However, well within 

L R L Ft L R 

the Coulomb diamond we have Fjv<~jv ^ ^dk' ^^dk ^^d'd' ^^'l perturbation theory can be used. This is the regime 
of interest here. 



IV. INELASTIC CO-TUNNELING 



If the bias voltage can supply_the necessary excitation energy, > A, then the excited state can be populated 
by inelastic co-tunneling events. □ In such a coherent process, an electron is transferred from the left to the right 
lead, thereby loosing an amount A of its excess energy, which is supplied by the transport voltage; this energy is 
transferred to the dot, which is thereby excited. The excited state can relax to the N electron ground state by two 
sorts of processes depending on A. If A is sufficiently small then there is an energy barrier for injecting or emitting 



electrons [i.e. conditions ( [7 a)) and (7b) are not satisfied]. Thus the excited state is also Coulomb blockaded. The dot 
can relax by three different inelastic co-tunneling processes. An electron can be transferred from the left to the right 
lead again, but now gain an energy A, thereby relaxing the dot. Similarly, an electron injected from either lead can be 
excited above the Fermi level of that same lead. The dot relaxes without a net charge transfer. This is the situation in 
the region labeled (I) in Fig. ^ For the current the elastic co-tunneling through the excited state is also important. If 
A is sufficiently large, however, the excited state can decay to the ground state by two sequential tunneling processes. 
There are three possibilities. Either the dot can relax to the iV — 1 electron ground state and then return to the N 
electron ground state [region (II) in Fig. ; or the dot can be excited to the N + 1 electron ground state from which 
it returns to the N electron ground state [region (IF)]; or both processes can occur [region (III)]. For positive bias 



voltage these regions are delimited by the bias voltage threshold lines [Eqs. (^) and (^] as follows 

(I) A <V< v^'^-W^^'"^' A < 1^,, 

(II) A, <V< Vr-\Vr'''' A < ^^ax 

(IF) A,Vf +1^' <V < vl!+^^,V^'^-^ A < iCa, 



.. .. 2 max 

(III) A,y/^-i,^/+i^' <v < v^'^-\v^+^'' A < C 
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Depending on the excitation energy A = 5 — J relative to the maximum bias voltage VJ^^ = S + U some of these 
regions do not exist. The different limiting cases are depicted in Fig. ||. When the excitation energy is sufficiently 
small, A < Vjj^x/S, all these regions exist. From this we find that the Coulomb and/or exchange interaction should 
large relative to the level spacing: S < + ^J. Region (I) ceases to exist for larger values of A (or i5) while regions 
(11) and (IF) continue to exists as long as A < i.e., S <U + 2J. Region (III) always exists since A < V^^^ is 

satisfied for U + J > 0. The situation is clearly more complicated than at low voltage y < A. The transport in each 
region needs to be considered separately. Below we first consider the simplest region (I), where only the ground and 
excited state participate in the transport. Then we extend this analysis to the other regions. 



A. Region I - relaxation by co-tunneling 



This region exists only when the excited state lies sufficiently close to the ground state. There are four states 
involved in the transport, of which three are degenerate. In order to calculate the current we need the nonequilibrium 
occupation of the states in the dot. We assume that at the bias voltage where inelastic co-tunneling processes become 
important, V ^ A, the lifetime of the states in the dots is much smaller than the typical time scale for building up 
Kondo-type correlations between the leads and dots system. In terms of energies we thus assume Tr- <^ A < ^ F 
where Tk is the energy scale for such correlations. In this limit we can use the following master equationalj to describe 
the occupation of the ground state pnq o a-^d of the three triplet excited states p^^ ,m = 0, ±1: 



dtp No 



= 1,2,3 



-N' 



J2 p^i. 

1=0, ±1 



3F 



N'^NPNa^o^ 



(14a) 



(^tPN[ ^ - ^N'^NPNo,o + ^r^'^AT/ ^ PN^, 



m=±l 



i=l,2,3 



PN' 



2^N'^N'^ ^ N*-N' 

1=1,2,3 



(14b) 



(14c) 



where m = ±1 and pNa o + Sm=o ±i PN'i ~ ^- current reads 



/ — {2Tn^n + SFjv'^at) pno,o + 



(1) 



-N' 



.(1) 



m=0,±l 



(15) 



The transitions between the states are schematically indicated in Fig. ^. The rates are calculated using the Golden 
Rule as was done for Eq. (^2|) . In Table | we have listed the final state of the leads, the matrix element and the 
transition rate for each transition. The explicit expressions for the matrix elements are given in Table ||. Let us now 
calculate these rates and discuss the transitions that they describe starting with the inelastic co-tunneling processes. 

Each excited state of the triplet is populated by inelastic co-tunneling transition from the ground state with the 
same rate Fjv^at. An electron with energy ei and spin s — ±1/2 is first injected into the dot from lead L, bringing 
the dot in the intermediate state |iV -I- l)i/2,si and is subsequently emitted into the lead R at an energy Sr = Si — A. 
The other sequence is also possible, in which case the intermediate state of the dot is \N — l)i/2,-is- The sum of the 
amplitudes for these processes yields the matrix elem ent M js'^n, listed in Table The inelastic co-tunneling rate 
for exciting the dot is then reads [cf. Eqs. (13a) and ( 13b )] 



= 2^ ^ \Mm'^n\^ f [ei - ml) [1 - / (sr ~ Pr)] S {ei - - A) 



Lr 



2tt 



^d'd'^dd [r {Enn-1, V-A)+t {En+in - A,V - A) 



-2t' (Enn-i, En+in ~A,V- A)] 



(16a) 



(16b) 



Conversely, each state of the triplet can relax to the singlet ground state by inelastic co-tunneling. The rate for 
these processes must be calculated separately because they result in a different final state of the leads. The matrix 
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element M^^^, in Table || describes the transfer of an electron from the left to the right lead starting in any of the 
three triplet sublevels of the excited state. The energy of the emitted electron is increased by the excitation energy 
= £/ + A. The inelastic co-tunneling rate for this relaxation process is [cf. Eq. (^)] 



/ (ej - ^i) [1 - / (e, - Aifl)] (5 (ej + A - e,) (17a) 



= ^^i'A V {En'n^i, V + A)+t {En+in' +A,V + A) 

+2t' {En'n-i,En+in' +A,V + A)]. (17b) 

This rate explicitly depends on the voltage because the voltage window (both fj.L and fj,R appear in the Fermi 
functions). The second process involves the injection of an electron at energy ei from the left lead and the emission 

(2) 

of an electron to an empty state in the same lead at a higher energy ei' = e/ + A. The matrix element Mj^^^^, is thu s 
found by replacing r I ' , cf. Tabic ||. The tunneling rate is similarly found by replacing r, R ^ I' , L in Eq. (17a) 
and from the result (17b) we obtain: 

p(2) _ £iiT.(l) fy _ Q^ 
^ dd 

The third process involves the injection and emission of an electron to and from the right lead: replacing L —i- R we 
find: 

p(3) _ ^d'd' p(l) (y _ Q^ 
d'd' 

Let us now consider the elastic co-tunneling processes in which electrons are transferred from the left to the right 
lead. These involve four spin-degenerate intermediate states of the dot. In addition to the four ground states (Q), 
there are four other intermediate excited states 



I (N - l)'l/2, + l/2> - «i'T|0>, I (N + l)'l/2, + l/2> 

I (N - - at,jO), I {N + l);/2,_i/2) 



^ ^d'1^d'l^d] \ 



The addition energies for the transitions from the N electron excited state to these states are 

En'N-V — EnN-1 + a — (5 = E]yN_i + J, 

En+1'N' = En+in + S — a = En+in — J- 

These energies lie closer to the Fermi energies than Enn^i and En-^-in. Therefore the contributions of these in- 
termediate states to elastic co-tunneling amplitudes are larger than those in the amplitude Mn^n and they cannot 
be disregarded. Three types of elastic co-tunneling processes must be considered separately. Two of these processes 
do not involve a flip of the spin of the transmitted electron. The two different final states of the leads as indicated 
in Table |. The first type of process occurs with the same rate ^rj^)^^, for all states of the triplet. The second 
type is only possible for the two states with nonzero spin projection and occurs with a rate ^rj^)^^,. The rates are 
calculated from the matrix elements M^^/^jv' Table || and give [compare this with Eq. (|lS 

J- N'^N' — ^TT > \1U 



Lr 

1 '^L T^R 



f {ei - ^1L) [1 - / {sr - fin)] S [e, - Sr) (20a) 

o {^dd^M [r iEN'N-i,V) + T {En+in',V) ± 2t' {En' n-i, En+in' ,V)] 

+^d'd'^d'd' [t {En'n-v.V) + T {En+vn',V) ± 2t' {En'n-v ,En+vn' .V)] 
+2 ReF^^/Pf,^ [r' {En'n-i,En'n~i': V) ± r' (i?jv'jv-i, £^Af+i'A" , V) 

+t' {En+in',En+i'N',V) ± t' {En'n-i'.En+in',V)]} . (20b) 

These processes do not induce transitions between the states [they don not contribute to Eqs. (|4|)], but they do 
contribute to the transfer of electrons [the rate appears in Eq. (Esl)]. The third type of elastic co-tunneling process, 
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involves a flip of the spin of the transmitted electron and induces transitions between the triplet states states. The 



triplet state |A^{_o) can decay to both |-/V{_i) and \N{ _i) with a rate equal to ^T^)^^,. The latter two states can 



return to the former state with the same transition rate. 

Because there is no Zeeman energy or spin-orbit coupling in our model, the transition rates between the states 
of the triplet balance each other, cf. Fig. ^. Therefore we can introduce the total occupations of the excited state 
PN' = J2m=o ±1 Pn[ summed over the spin projection and omit the spin indices pN = PNq o- By adding the equations 
for the triplet states we obtain a closed set of master equations for the total occupations: 

dtPN = ^ ^'^m^m'Pn' ~ STn'^nPn, (21a) 

i=l,2,3 

dtPN' = S^N'^nPn - ^ ^n^^n^Pn', (21b) 

1=1,2,3 

with pM + PN' = 1. The elastic co-tunneling transitions between the triplet states have precisely cancelled. However, 
these processes do contribute to the current (p5|). The fact that the triplet state with zero spin projection transfers 

twice as much charge per unit time (2 x ^T^'^^,) as the other states (|r^)^^,) is precisely balanced by the spin-flip 

processes that can only occur for states with nonzero spin projection (irj^)^jY,). The elastic co-tunneling contribution 
of each state of the triplet is thus proportional to the same rate (ttq) and we can also express the current in total 
occupations 



/ = {2Tn^n + arjv'^w) PN + ( ^r^^^, + r^L^v + ^nIn' ] pn'- (22) 



Note that the inelastic co-tunneling relaxation processes with rates rj^^jy, do not appear here since they do not 
involve charge transfer. In the stationary limit limt^oo dtPN' , dtPN — equations ( pT| ) are readily solved. The current 
can be written as 

/'i-r(i) I r(2) _ or ^ .l or^^) -l r(*) 
I — ^-^ "STn'^n + '^^N^N- (23) 

Si=l,2,3 ^N'^N + Sr^v'^AT 

For < A the excitation rate Fat'^at vanishes and we recover the elastic co-tunneling current 111]. In Fig. |4| we have 
plotted the current, differential conductance and the co-tunncling rates as a function of bias voltages for the case 
A < V^^Ji. In this hmit region (I) covers most of the Coulomb diamond for ^ > A. Apart from a clear kink in the 
current and a step in the differential conductance, a sharp peak is present at the step edge. This is clearly due to a 
change of the population of the excited state. 



B. Regions (II), (II') and (III) - relaxation by sequential tunneling 

The remain three regions where sequential tunneling transitions from the excited state to — 1 and/or + 1 
electron ground state are possible are now easily included. The elastic co-tunneling processes involving the excited 
state must now be disregarded ( and their contributions to the current). The rates for these processes diverge at 
the boundaries of regime (I) signalling that perturbation theory breaks down. Instead we write a full set of master 
equations similar to (|l^) without the transitions between the triplet states and now including the occupations of the 
spin sublevels of the A^ — 1 and/or A^ + 1 electron ground states. Proceeding as before we arrive at the following 
closed set of master equations for the total occupations, including pn±i = X]m=±i/2 P{n±i)^^2 ' 

dtPN~i = a^Tf,^,pN' - (a-T^j^ + a+r^,^,) pjy-i, (24a) 

dtPN+i = a+TddPN' - {a-^d'd' + o^+'^d'd') Pn+i, (24b) 

dtPN = +o:-TjdPN-i + a+T^i^ipN+i - 3Tn'*~nPn, (24c) 

dtPN' = ct-^d'd'PN-i + a+Tf^pN+i + 3Tn'^nPn 

- i+o^-'^d'd'PN' + «+rL-OA") PN'. (24d) 

The current reads 
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/ = (2r^^Ar + 3Tn'^n) Pn + a.T^,^,pN' + «+ {T^d + ^d'd') Pn+i- (25) 

These expression covers regions (II) [q!_ — 1,q;+ — 0], (IF) [a_ — 0, = 1] and (III) [a± — 1]. The solution for all 
cases can be written in the form 

A + 6K,l N'^N 

The coefficients appearing here are constants that characterize the differences between the three regions: 

^ dd I „ dd 

K ^ a-IrfM' T^r , r.L ^^+'-d'd' 



pL I -pL ' + d'd' pfl, I pfl : 
^ d'd' dd d'd' 

TV-' „ pi?, ^ d'd' I „ pL ^ d'd' 



d'd' pL I pi ^ + rf'd' pi? I pfl ' 
^ dd d'd' dd d'd' 

pi?. pi 

/V = a- , +a+ „ + 1. 

pL I pL T pfl I -pj^ 

^ dd ^ d'd' dd ^ d'd' 

where K + K' — a_r^^, + a+F^,^, . In Fig. ^ we have plotted the co-tunneling rates, the occupations of the states 
and the differential conductance as a function of the bias voltages for T^^ — F^,^, = F, a = L, i? for a typical case 
where < A. In this limit only region (III) exists. The qualitative result is the same as in region (I). The 

essential condition for the peak in the conductance is that the total rate for charge transport through the excited 
state is bigger than the rate F^v't-jv at which this state is populated. 



V. CONCLUSIONS 



We have analyzed the finite bias transport within a Coulomb blockaded region of a quantum dot with an even number 
of electrons. The differential conductance through the dot was found to exhibit a step within the Coulomb diamond at 
a bias voltage equal to the excitation energy. At the step edge an additional peak appears that is caused by a change 
in the nonequilibrium populations of the dot states due the voltage dependence of the co-tunneling rates. This may 
explain the recent observations in the transport spectroscopy of serai-conductor quantum dots by De Franceschi et. al 
and of carbon nanotubes by Nygard et. al. As pointed out in Ref.El the sharp feature in the differential conductance 
can be used to do a spectroscopy of excited states with increased resolution. We have considered two well-separated 
regimes where the excited state is depopulated by inelastic co-tunneling and sequential tunneling, respectively. It is 
expected that in the crossover regime the results do not change qualitatively since the inelastic co-tunneling processes 
that excite the dot form the bottleneck for the transport. The rate for these processes as calculated here is correct 
also in the crossover regime. The calculations presented here can be improved by including (additive) contributions of 
more virtual states of the dot (more than two orbitals) to the transition amplitudes, which determine the co-tunneling 
rates. Also, the situation where a second iV-electron excited state lies the Coulomb diamond can be treated in the 
same way as single excited state. The occupations of all transport states can be calculated from an extended set of 
master equations. One must now also include co-tunneling transitions between the two excited states. The rates that 
enter in these equations have the same form as those calculated here. Furthermore, the peak at the on^et of inelastic 
co-tunneling is broadened by thermal fluctuations and level broadening effects, as pointed out in Ref.cl. The effect of 
a electron finite temperature can be included in our calculations. Level broadening effects require a nonperturbative 
scheme, which is beyond the scope of this paper. It is clear, however, that near the onset of inelastic co-tunneling, 
the broadening of the peak, which is determined the decay rate of the excited state (either by sequential or inelastic 
co-tunneling), is much larger than the rate for exciting the dot. 
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TABLE I. Transitions from initial state |i)a[|.ajj) to final state with the second order matrix elements and corre- 

sponding rates. The explicit expressions for the matrix elements are given in Table |l| The rates are given in the text. Here I 
and r' label occupied states in the left and right lead, respectively; I' and r label unoccupied states in the left and right lead, 
respectively. 



type 
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f 






F 






elastic 


No,o 


No,o 




1 

V2 




V2Mn<-n 


2r]v^]v 


inelastic 


No,o 


Ko 




1 






Fat'^-jv 




No,o 
No,o 


K+i 

K-i 










elastic 


NU 




= ±1,0 


1 
V2 






ip(l) 

2-"- JV'^JV' 








= ±1 (!) 
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Ip(^) 

2^ N'<—N' 




Ko 
Ko 
K+i 


K + 1 

K-1 
Ko 










2-"- JV'<-JV' 


inelastic 


Ko 

K+i 

K-i 


No,o 
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No.o 




1 
V2 


(«iVrT + «?T«n) 

"it rT 
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Ko 
K+ 
K- 



Ko 
K+i 

K-i 



No,o 
No,o 
No.o 



75 K^I'T 



No,o 
No.o 



-T= (a^ , I a\ + a^-^a^ , ) 

„t „t 



M 



M 
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TABLE II. Second order matrix elements for the transitions listed in Table | . 



M^^j, = ^'f" + ^ elastic 

£r ~ iijvJV-1 iijv+iJV — £i 
j^(l,2) _ tldtrd ^id'^rd' 



£r — £r ^ -BjV'JV-1' 



En+in' ^ £; Ej^^irj^i — El 



Mj,,^^ = ^'"'^'-^ + ^ inelastic 

Er — EnN-1 En + 1N — £l 
^(1) _ ild^rd' I tld'trd' 



Er — SjV'JV-l -Ejv-i-ijv' — £; 



,^(2) _ ti^tl'd' ttdtl'd' 

£l' ~ -C'JV'JV-1 -C-JV-j^liV' — 
1^(3) _ tr'd^rd' iy'd^rd' 



E 



JV+IJV' 



FIG. 1. Stabilitity diagram of a quantum dot for the symmetric case Cl/Cq ~ CrJCg 3> 1. The skewed lines delimit the 
regions where sequential tunneling is possible through the ground states [thick; Eqs. (^)] and the excited state [thin; Eqs. (^)]. 
The horizontal thick lines mark the onset of inelastic co-tunneling. 



FIG. 2. Inelastic co-tunneling regions for the different cases of excitation energy relative to the maximal bias voltage. 



FIG. 3. Transitions between the spin-singlet ground state and the spin-triplet excited state. Elastic co-tunneling processes 
which do not change the state (charge and spin projection) are indicated by a circle, beginning and ending at the same state. 



FIG. 4. Results for region (I) as a function of bias voltage in the center of the Coulomb diamond; A/V^^^ = 0.1, 5/U — 5, 
U — 1, T/U — 0.05, Cl/Cg ~ Cr/Cg = 0.1. Left: Co-tunneling rates relative to the sequential tunneling ratee, i.e., in units 
of r. Elastic: (1) for the A'^ electron ground state Fjv^jv (2) for the A'^ electron excited state F^^^j^, = Tiv'^jv'- Inelastic: 
(3) for exciting the dot Fjv^jv (4) for relaxing the dot with charge transfer Fj^^j^, and (5) without charge transfer F^^'j^, . 
Middle: Occupations of (1) the A^ electron ground state (2), the A' electron excited state and (3) the A' — 1 and A' -f 1 electron 
ground states (coincide). Right: Differential conductance calculated from expression (p3|) for the current. 



FIG. 5. Results in region (III) as a function of bias voltage in the center of the Coulomb diamond: A/V^^^ — 0.5, S/U — 5 
,U = l,F/(7 — 0.04, Cl/Cg ~ Cr/Cg = 0.1. Left: Co-tunneling rates relative to the sequential tunneling rate, i.e., in units of 
F. Elastic: (1) for the TV elctron ground state Fjv^jv. Inelastic: (2) for exciting the dot Fjv'^jv- Middle: Occupations of (1) 
the A^ electron ground state (2) the A electron excited state and (3) the A — 1 and A -|- 1 electron ground states (coincide). 
Right: Differential conductance calculated from expression (E^) for the current. 



11 



Co-tunneling rates Occupations Differential conductance 




V/A V/A V/A 
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